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Ab initio simulations are combined with in situ infrared spectroscopy to unveil the molecular trans- 
port of H2, CO2, and H2O in the metal organic framework MOF-74-Mg. Our study uncovers — at 
the atomistic level — the major factors governing the transport mechanism of these small molecules. 
In particular, we identify four key diffusion mechanisms and calculate the corresponding diffusion 
barriers, which are nicely confirmed by time-resolved infrared experiments. We also answer a long- 
standing question about the existence of secondary adsorption sites for the guest molecules, and we 
show how those sites affect the macroscopic diffusion properties. Our findings are important to gain 
a fundamental understanding of the diffusion processes in these nano-porous materials, with direct 
implications for the usability of MOFs in gas sequestration and storage applications. 



PACS numbers: 68.43.Bc, 68.43.Mn, 84.60.-h, 84.60.Ve 

Distressing scenarios of fossil-fuel shortages and green- 
house effects have intensified scientific efforts in seeking 
innovative materials for energy storage and capture of 
greenhouse gases. Due to their vast range of tunable 
properties, metal organic framework (MOF) materials 
are outstanding candidates to address both issues. They 
can have large hydrogen-storage capacities, surpassing 
the US DOE storage targets [l--3 j , and have even been de- 
clared DOE's top priority at its 2012 Hydrogen and Fuel 
Cells Program Annual Review Meeting [3] . On the other 
hand, due to their large CO2 uptakes |B], MOFs are 
also interesting for carbon-capture applications. These 
qualities, combined with their low production cost, have 
made MOFs the target of many studies, focusing mostly 
on their adsorption properties [3] [3] [THS] ■ However, their 
performance for practical storage and capture applica- 
tions also critically depends on the diffusion and interac- 
tion of gases (e.g. with water) in the MOF environment, 
which are currently poorly understood [6] [TO] [UJ . 

To address this issue and elucidate the diffusion process 
of small molecules in the nano-pores of MOFs, we use a 
combination of ab initio simulations and in situ infrared 
(1R) spectroscopy. Specifically, we study the diffusion 
of H 2 , C0 2 , and H 2 in MOF-74-Mg and investigate 
the role of water on diffusion of small molecules. We 
focus on this particular MOF, as it has attracted a lot 
of attention due to its enhanced reactivity with small 
molecules, caused by the exposure of open metal sites in 
its nano-pores. Britt et al. [5] have shown that MOF- 74- 
Mg is extremely efficient in capturing CO2, compared to 
iso-structural MOFs with other metal sites such as Zn, 
Mn, Fe, Co, Ni, and Cu. 

To model the molecular diffusion in the MOF struc- 
ture, we use climbing-image nudged elastic band |12j 
(NEB) simulations coupled with density functional the- 
ory (DFT), as implemented in QuantumEspresso [13]. 
To correctly capture the weak van der Waals forces — 



which are critical in this study — we use the truly non- 
local functional vdW-DF [14TU6| . which has already 
been successfully applied to a number of related stud- 
ies [P7Hz!T] . We use ultrasoft pseudopotentials with wave- 
function and density cutoffs of 480 eV and 4800 eV, re- 
spectively. Tests show that T-point calculations are suffi- 
cient and yield total energies converged to within 5 meV 
with respect to denser fc-point meshes; however, en- 
ergy differences — important for our diffusion barriers — 
are converged to within less than 1 meV. The self- 
consistency tolerance was set to 1.4 x 10 -10 eV and dur- 
ing optimizations the total forces were relaxed to less 
than 2.6 x 10 -4 eV/A. We started from the experimental 
MOF-74-Mg structure, optimizing the internal parame- 
ters and keeping the lattice constants fixed to the hexag- 
onal structure with space group R3 and a = 25.8810 A 
and c = 6.8789 A [52]. 

We begin by revisiting the binding characteristics of 
H2, CO2, and H2O in the MOF framework, which 
serves here as background for explaining the molecu- 
lar transport — recent theoretical |23H25] and experimen- 
tal [2~4Tf2"8] work has fully explored and clarified both ori- 
gin and nature of the "static" adsorption interactions of 
H 2 and CO2 in MOF-74. Our calculated adsorption ener- 
gies AE are given in Table|l]for two different loadings: (i) 
low-loading, one molecule per cell, and (ii) high-loading, 
six molecules per cell completely saturating all available 
metal sites (see Fig. [T^, and Fig. [T^,' for a graphical rep- 
resentation of these loadings). From the computed vi- 
brational frequencies (see below) we obtain the thermal 
and zero-point energy (ZPE) corrections to the adsorp- 
tion energy, allowing for a more accurate comparison to 
measured adsorption heats. But, we find that neither of 
these corrections has a large effect in our case. Over- 
all, we find very good agreement with the experimental 
adsorption energies of -0.11 ± 0.003 eV for H 2 [27] and 
-0.49 ± 0.010 eV for C0 2 [21], attesting to the impor- 
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TABLE I. Adsorption energy AE of molecules in MOF-74-Mg 
in eV, and adsorption energies corrected for the zero-point 
energy (A_Ezpe) and thermal contribution (A//298 at 298 K). 



Molecule 


Loading 


AE 


ASzpe 




H 2 


1 


-0.15 


-0.15 


-0.15 












6 


-0.16 


-0.16 


-0.16 


C0 2 


1 


-0.50 


-0.49 


-0.50 




6 


-0.50 


-0.49 


-0.50 


H 2 


1 


-0.79 


-0.76 


-0.76 












6 


-0.76 


-0.73 


-0.73 



tance of correctly including van der Waals interactions in 
these simulations. The 0.03 eV decrease of AE for H2O 
in the high-loading situation is linked to the intermolecu- 
lar repulsions between hydrogen atoms of H2O adsorbed 
on adjacent Mg sites. This is also evident from the inter- 
molecular distance dHHO - Mg-MOF, which increases from 
2.197 A for the low- loading case to 2.232 A for the high- 
loading case. As a result of the large dipole moment 
of water (1.9 Debye), it binds very strongly [2J. Thus, 
water is thermodynamically much more likely to occupy 
metal sites, significantly reducing the MOF's adsorption 
capabilities towards H2 and CO2, with tremendous im- 
portance for storage and capture applications. 

While the uncoordinated Mg sites are the most attrac- 
tive sites, secondary binding sites also exist, which only 
become occupied at lower temperatures (or high pres- 
sures). Recent neutron diffraction experiments isolated 
such a secondary binding site for CO2 in MOF-74-Mg 
[28] . where the CO2 binds to the carboxylate group of 
the linker. The experimental geometry of this site is well 
reproduced by our vdW-DF calculations [25 . For very 
high loading, i.e. 12 CO2 molecules per unit cell (occupy- 
ing the primary and secondary sites) the inter-molecular 
repulsions decrease the average AE to -0.48 eV. 

Vibrational frequencies are a prerequisite to estimate 
the pre-exponential factor in the Arrhenius transport 
equation. We report here the calculated change in 
IR frequency between the gas-phase molecules and the 
molecules adsorbed in the MOF. In the mono-adsorbed 
cases we find A^h 2 = -30 cm -1 , Avco 2 — "13 cm -1 , 
and Az^h 2 o = -103 cm -1 , in excellent agreement with 
our experimental observations of -36 cm -1 , -8 cm" 1 and 
-99 cm . Experimentally, there is a small difference 
in frequency change between low and high loading, re- 
sulting in a red-shift of -3 cm" 1 and -15 cm -1 for the 
asymmetric stretch mode of CO2 and H2O (see supple- 
mentary material). Computationally, we find Av>co 2 — 
-1 cm -1 and Afn 2 o = _ 18 cm" 1 , in excellent agreement 
with experiment. 

We now move our discussion to the key results of this 
letter: The molecular transport in MOF-74-Mg. To 
study this aspect, we consider the four diffusion mech- 
anisms depicted in Fig. fl] In mechanism a) , a molecule 



TABLE II. Diffusion barriers AE D (in eV) of small molecules 
in MOF-74-Mg. ZPE corrected barriers AE% PE are given in 
brackets. 







diffusion mechanism 






a) 


b) 


c) 


d) 


H 2 


0.07 (0.07) 


0.15 (0.15) 


0.005 (0.004) 


0.16 (0.16) 


C0 2 


0.27 (0.26) 


0.29 (0.28) 


0.04 (0.03) 


0.45 (0.44) 


H 2 


0.67 (0.64) 


0.61 (0.58) 


0.06 (0.04) 


0.73 (0.70) 



hops circularly from one Mg 2+ site to its adjacent one. 
Note that this mechanism is not responsible for molec- 
ular transport into the MOF, but nevertheless it is an 
important process for redistributing the load. In mech- 
anism b), one molecule hops longitudinally (along the 
c-axis) from one Mg 2+ to the equivalent one in the next 
unit cell. In mechanism c), one molecule moves longi- 
tudinally through MOF- 74 fully loaded with the same 
type of molecule. And finally, in mechanism d), one 
molecule pre-adsorbed on a Mg 2+ site moves through a 
barrier made by six molecules and then binds again at 
the equivalent site two unit cells further down. We con- 
sider those the fundamental diffusion mechanisms that 
control the macroscopic molecular transport in MOF- 74- 
Mg. Other paths by good approximation are superpo- 
sitions of the ones discussed here. Mechanism c) sim- 
ulates real diffusion and is responsible for penetration 
deep into the MOF, once the surface is fully saturated. 
Mechanism d) simulates the kinetic barrier that a pre- 
adsorbed molecule must overcome when the adjacent 
available metals sites have already been saturated, i.e. 
obstructing the flux of an incoming molecule. Note that 
co-diffusion cases are not considered here. For the study 
of mechanism a) we used the hexagonal primitive cell 
with 54 atoms. For mechanism b), c), and d) — requiring 
a longitudinal displacement of the molecule along the c- 
axis — we employed a supercell containing 108 atoms, ex- 
panded along the c-axis, making these calculations very 
challenging. 

Figure [2j consistently labeled with Fig. [lj shows the 
calculated diffusion barriers for H2, CO2, and H2O. The 
barrier heights of all curves, i.e. their maxima, will be 
referred to as AE D and are collected in Table Hll Not 
surprisingly, looking at the circular hopping mechanism 
a), we see that the molecules experience diffusion barri- 
ers similar in magnitudes to the adsorption energies at 
the metal site. For both H2 and H2O the diffusion bar- 
riers look alike and are found very symmetric. This is 
not the case for CO2, which shows a second minimum 
located at 58% (see Fig. [2^). The local geometry of 
this minimum is that of the secondary adsorption site 
detected by neutron diffraction in Ref. EH Figure [2^l) 
shows that the secondary adsorption site becomes only 
activated when all the six Mg 2+ sites are fully occupied, 
as the local minimum (5 meV) is too shallow to trap 
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FIG. 1. Graphical representation of the diffusion mecha- 
nisms considered in this study, shown for the case of CO2. a) 
and a') are views directly along the c-axis of the hexagonal 
MOF-74-Mg cell, where one (low loading) and six CO2 (high 
loading) are adsorbed, b), c), and d) are views perpendicular 
to the c-axis. Dashed lines indicate the diffusion paths. 



molecules even at very low temperatures. In other words, 
the secondary binding site only gets occupied under very 
high loading/pressure. Barriers calculated for mechanism 
b) are similar in magnitude to those of mechanism a), 
however, they are slightly more asymmetric. We con- 
sidered mechanism b) also for situations of low-loading 
(not shown in Fig. |2|, i.e. only one molecule occupies 
the unit cell and hops along the c-axis. In this case, 
all barriers remain similar with the exception of that for 
water, which increases by approximately 10%, suggesting 
that the crowded environment helps the molecule extrac- 
tion via an expulsion mechanism. Not coincidentally, the 
maximum of the H2O barrier of Fig. ^p) corresponds to 
the formation-breaking of a lateral hydrogen-bond with 
an oxygen nearby the Mg sites (OOH- • -O-MOF), which 
decreases its strength passing from low (dooH - o-MOF = 
2.077 A) to high loadings (rfoOH— O-MOF = 2.138 A). It 
follows that H2O experiences the largest diffusion barri- 
ers and has the highest adsorption affinity to metal sites 
(see Table |l|. 

Although not immediately obvious from Fig. [2]d), for 
H2 there is a small secondary minimum near the 50% 



FIG. 2. Diffusion profiles (in eV) for the diffusion processes 
of H 2 , C0 2 , and H 2 in MOF-74-Mg according to the mech- 
anisms in Fig[l] 



mark, which is off the optimized pathway and thus does 
not show up in the NEB calculation. With the help of 
additional optimizations near that point, we confirm the 
presence of a secondary binding site with an adsorption 
energy of AE = —0.11 eV, which is 0.04 eV above the 
primary binding site (see Table E|. The binding pocket is 
relatively shallow with a depth of approximately 2 meV 
and the H 2 is clearly "docked" on the "oxygen triangle" 
formed by the linkers coordinating the Mg sites. Evi- 
dence of this binding site was also found in MOF-74-Zn 
both theoretically |17j and with neutron diffraction ex- 
periments [29 . Our calculations thus confirm and quan- 
tify the schematic kinetic model (based on experimental 
evidences) proposed in Fig. 12 of Ref. |Hj, where the 
gas-phase H2 molecules first condense on the metastable 
second binding site, followed by a gradual migration to 
the more stable Mg site. Also note that the asymmetric 
profile for this mechanism suggests that the release of a 
H2 molecule from the Mg site is a direct process with- 
out intermediate steps. From this analysis we learn that 
at high pressure and temperature the existence of sec- 
ondary binding sites assists in the quick redistribution of 
the molecular load, as molecules entrapped in these sites 
need less energy to escape from their minima. These sites 
thus play a crucial role in the load redistribution, affect- 
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ing the final diffusion properties of incoming molecules. 
For example, when water is adsorbed at the primary and 
secondary sites, we can imagine scenarios of high hydro- 
gen bonding reticulation, severely limiting the molecular 
transport of H2 and CO2. 

Next, we focus on mechanism c) and d), which aim 
to simulate the molecular flow in the MOF. In mecha- 
nism c), we calculate the barrier for diffusion of a sin- 
gle molecule through a MOF fully loaded with the same 
type of molecule. Again, water suffers the largest imped- 
iment to flow (0.06 eV), while H2 for this mechanism is 
almost barrier-free (5 meV). Just recently, Bao et al. [5] 
measured AE D for CO2 and found 0.03 eV, based on 
experimental CO2 adsorption isotherms. This is in per- 
fect agreement with our computed values of 0.04 eV and 
0.03 eV (ZPE corrected). Mechanism d) considers a situ- 
ation where an adsorbed molecule tries to diffuse further 
into the MOF, hampered by the presence of a molecular 
barrier due to other pre-adsorbed molecules. In this case, 
the strong electrostatic interaction established by water 
with the Mg sites imposes a larger diffusion barrier for 
H2O when compared to CO2 and H2. 

We now correlate our diffusion mechanisms at the 
atomistic level discussed above with our in situ time- 
resolved IR diffusion experiments, which provide infor- 
mation about the macroscopic diffusion of molecules 
within the MOF scaffold. To this end, we performed IR 
measurements with diffuse reflectance geometry using a 
liquid-N 2 cooled Indium Antimonide detector. Approxi- 
mately 12 mg of sample was mixed with 0.02 g of KBr and 
inserted into an atmospheric cell from Thermo scientific. 
The sample was heated in situ at 450 K in vacuum for 
complete de-solvation. The pressure was kept constant 
and changes were monitored as a function of time, while 
temperature was set to 298 K. IR spectra were recorded 
every 10 seconds. 

By measuring the asymmetric stretch modes of both 
C0 2 and H 2 with time, we find that typical MOF load- 
ing occurs in two consecutive steps. At the beginning, 
the guest molecules condense in nano-pores on the sur- 
face of the MOF, leading to very high- loading in those 
nano-pores. From our experiments we conclude that this 
corresponds to approximately 9 molecules per unit cell, 
effectively clogging the pathways and limiting diffusion 
of new approaching molecules [55] . This results in a red- 
shift (see supplementary materials), explained by increas- 
ing lateral interactions between adsorbed molecules |25j . 
Our own calculations of the frequency shift between low 
and high loading from above confirm that this shift cor- 
responds to high loading. As time progresses, molecules 
diffuse from highly-loaded pores to low-loading pores 
deeper inside the bulk of the MOF. As a result, the 
asymmetric mode returns (with a blue-shift, see supple- 
mentary material) to its original location. This peculiar 
fluctuation in frequency-shifts is a clear signal of a two- 
states molecular transport mechanism. In particular, the 



time it takes to return to the original shift — i.e. to get 
from high to low loading — is characteristic for the guest 
molecules and provides insight into the diffusion process. 

We have performed the above described diffusion ex- 
periment for CO2 and H2O, and find that it takes ap- 
proximately 2 hour for water to diffuse to the low-loading 
situation, while it takes only 22 minutes for CO2 (due to 
the difficulties of the experiments, we report these times 
with a large error bar of approximately 25%). More 
importantly, the ratio of those two times is 5.45. The 
described experimental situation clearly corresponds to 
diffusion through mechanism c) and we can calculate the 
corresponding ratio using the Arrhenius equation. To 
this end, we estimated the pre-exponential diffusion fac- 
tor using the harmonic transition-state theory (see sup- 
plementary material) |30| . Using our barrier heights from 
Table |ll| for CO2 and H2O for mechanism c), we find for 
the same ratio a value of 5.43, in excellent agreement 
with experiment. This confirms that we have success- 
fully uncovered transport processes at the atomistic level, 
that govern the transport processes at the experimentally 
measured macroscopic level. Ideally, we would like to re- 
port other ratios, e.g. for H2/CO2, but the corresponding 
IR experiments involving H2 are very difficult and results 
will be published in a longer, forthcoming article. 

The fact that water diffuses so much slower than H2 
and CO2 is to be expected, but of paramount importance 
for practical applications of MOFs. While the presence of 
water does not hinder the initial diffusion of other gasses 
such as H 2 and CO2 in the non-contaminated MOF struc- 
ture, high pressure or temperature can affect the delicate 
thermodynamic equilibrium as incoming H2O molecules 
displace pre-adsorbed H2 and CO2 molecules. 

In summary, we demonstrate that state of the art ab 
initio NEB simulations, coupled with time-resolved in 
situ IR spectroscopy, provide a complete atomistic un- 
derstanding of the diffusion processes of three important 
molecules in MOF-74-Mg for the purpose of hydrogen 
storage and gas separation. Our theoretical atomistic 
model for the molecular transport explains experimen- 
tal IR macroscopic evidences. Furthermore, our calcula- 
tions clarify the two-state mechanism, observed exper- 
imentally, which controls the macroscopic diffusion of 
these molecules. While the present study only focuses 
on MOF-74-Mg, it is our belief that the same methodol- 
ogy can be successfully applied to unveil the molecular 
transport in other MOFs and nano-porous materials, pro- 
viding further insight into the important question of dif- 
fusion as well as a robust theoretical foundation to guide 
the interpretation of challenging diffusion experiments. 

This work was entirely supported by the Department 
of Energy Grant No. DE-FG02-08ER46491. 
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